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ABSTRACT 

If the density distribution p{r) of MACHOs is spherically symmetric with respect to 
the Galactic center, it is shown that the minimal total mass M^^ CHO of the MACHOs is 
1.7 x 10 10 Mqt^P where rj^-p is the optical depth (r LMC ) toward the Large Magellanic 
Cloud (LMC) in the unit of 2 x 10~ 7 . If p(r) is a decreasing function of r, it is proved that 
M^^ CHO is 5.6 x 10 10 Mqt^^. Several spherical and axially symmetric halo models of 
the Galaxy with a few free parameters are also considered. It is found that M^ CHO 
ranges from 5.6 x 10 10 Mqt l ^ to ~ 3 x 10 11 Mqt^ . For general case, the minimal 
column density S^cho of M ACHOs is obtained as E^acho = 25M pc- 2 r^c . if the 
clump of MACHOs exist only halfway between LMC and the sun, M^ CHO is 1.5 x 
1O 9 M . This shows that the total mass of MACHOs is smaller than 5 x 1O 1O M , i.e. ~ 
10% of the mass of the halo inside LMC, either if the density distribution of MACHOs 
is unusual or r LMC «2x 1(T 7 . 

(11 October 1996, Accepted for publication in Apj. Letters.) 
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1. Introduction 

Recent second year analysis of microlensing events 
toward the Large Magellanic Cloud (LMC) by the 
MACHO collaboration suggests that the optical depth 
t lmc is _ 2 x 10" 7 and the fraction / of MACHOs 
is ~ 0.4 with the typical mass ~ 0.34Mp) in the stan- 



dard sph erical flat rotation halo model ( Bennett ei 



al. 1996). The estimated mass of MACHOs is just the 
mass of red dwarfs. However the contribution of the 
halo red dwarfs to MACHO events should be small 



too low ( 


Bahcall et al. 1994 


praff & ] 


7 rccsc 1996b). 



Graff & Freese 1996a 



As for the white dwarf Galaxy halo, Chariot & Silk 
(1995) combined population synthesis models with 
constraints from deep galaxy surveys and showed that 
only a small fraction (< 10%) of the dark mass in 
the present-day galaxy halo could be in the form of 
white dwarf remnants (WDR) of intermediate-mass 
stars. Adams & Laughlin (1996) recently argued the 
implications of white dwarf Galactic halos. From the 



current limits on the densitv of red dwarfs ( 


Bahcall 


et al. 1994 


, Graff & Freese 1996a 


, Graff < 


ii Freese 



L996b ) and the galactic metallicity, the IMF must 
oe sharply peaked about a characteristic mass scale 
M c ~ 2.3Mq . They concluded that the mass frac- 
tion of WDRs in the halo is likely to be less than 25% 
since only a fraction of the initial mass of a star is 
incorporated into WDRs. 

The spatial density distribution function of MA- 
CHOs which caused microlensing events is not known 



of halo dark matter ( 


Paczvnski 1986 L Griest 1991 


Sackett & Gould 1992, 


Fricman & Scoccimarro 1994, 


sahu 1994 J Alcock et al. 1995, Gates, Gvuk & 


turner 199(17 turner, Gates, & Gvuk 1996, nivalis 


1996, 


Kan-ya, JNishi fc JNakamura 1996). Only the 



possible val ue of the optical depth r LMC ( Bennett 



et al. 1996) is known, the fraction of MACHOs in 
the halo depending on the spatial density distribution 
function of MACHOs. In this situation it is important 

to check the relation between the total mass M MACHO 
of MACHOs to the density distribution function of 

MACHOs for the given optical depth T LMC . The re- 
sults of such a study will be useful in the arguments 
on the fraction of MACHOs in the halo and what 
MACHOs are. 

In this paper we study jVf MACHO for various den- 
sity distribution functions of MACHOs and discuss 
the minimal total mass of the MACHO halo. Gates, 
Gyuk, and Turner (1996) and Turner, Gates, and 
Gyuk (1996) also discussed the mass of the MACHO 
halo for various galaxy models. While they separately 
added unidentified dark thick disk components, we 
discuss the total mass responsible for microlensing. 
In §2 we obtain the minimal total mass M^^ w of 
MACHOs for spherically symmetric density distribu- 



tions. In §3 we discuss M^ ACHO for various axially 
symmetric density distributions. §4 will be devoted 
to discussions. 

2. Spherically Symmetric Halo Models and 
the Minimal Total Mass of MACHOs 

We assume that the density distribution function 
p(r) of MACHOs is a function of the galactocentric 
radius r. The optical depth t lmc tow ard LMC is 



given by ( Paczyhski 1986, Griest 1991) 



LMC 



4ttG 



D. 



■)p(r)dx, (1) 



i?o - 2R j]x + x 2 , 



and 



(2) 
(3) 



T) = cos b cos I, 

where D s ,l,b and Rq are the distance to LMC (50kpc), 
the galactic longitude and latitude of LMC and the 
galactocentric radius of the sun (8.5kpc), respectively. 
In Eq. (1) we assumed the threshold «t=1 for sim- 
plicity. 

Equation (1) is rewritten as 



LMC 



/(*) 



and 



f(x)dm, 



dx 



dr 



dm = Anr p{r)—dx 
dx 



(4) 



(5) 



(6) 



For LMC, f(x) is infinite at x = x c = Rqi] = 0.153i?o 
so that the minimal total mass Af„ ACHO of MA- 
CHOs is zero for any given r LMC if MACHOs are 
distributed in an infinitesimally thin shell at r — r c = 

y/l — rfRo . However this is wrong. Since the angular 
size of LMC is - 10° x 10°, m ma cho is mm i m i ze d 
if MACHOs are distributed in a shell at r = r c with 
width d given by 



d = ^RoV - 227pc. 

It is easy to show that M^^p 110 is given by 

(r c + d f - rl 



(7) 



MACHO 



C 2 T LMC 



3G 2V 7 2^TW{R^-^L^ ) ' 



1.7 x lO lo M r^ 7 c , 



(8) 
(9) 

where is r LMC in the unit of 2 x 10~ 7 . This 

shows that in principle M MACHO can be only ~ 3% 



2 



of the total mass of the halo inside LMC. However 
the density distribution function of MACHOs in this 

case is very peculiar so that we calculate M MACHO for 

more realistic p(r) to know more realistic M^f ACH °. 

We consider two models; 

l)Polytropic Model 

p(r) is given by polytrope of index N and the radius 



Rv 



2) a Model 
p(r) is given by 



p(r) 



Po 



(1 + ^) 



(10) 



This model is similar to the beta model of the cluster 
of galaxies with core radius Ra- 
in figure 1 we show M MACHO as a function of 
R p for several polytropic indices N. M^A CHO ranges 
from 5.6x W w M q t^ for N=0 to 7.8X lO lo M r^ 
for N=3. For N=4 and 4.5, M MACHO is greater than 
1.0 x 10 11 Mqt^j and the minimum does not ex- 
ist for R p < D s . Under the assumption that p(r) 
is a decreasing function of r, it is shown in the Ap- 
pendix that M MACHO is minimized when p{r) is con- 
stant. Therefore M^ ACHO is 5.6 x lQ lo M rf^ if 
p(r) is a decreasing function. In figure 2 we show the 
total mass of MACHOs inside LMC in a models as 
a function of R a for several values of a. M^ ACHO 

ranges from 1.3 x 10 11 M@t^$ for a = 1.5 to 
8.7 x 1O 1o M tLmc for a = 6 For large ^ M^ ACHO 

does not change so much and it converges although 
the value of R a at the minimum increases . This 
behavior can be understood analytically using the 
asymptotic expression of gamma functions. 

3. Axially Symmetric Halo Models and the 
Minimal Total Mass of the MACHOs 

There are several suggestions that the Galactic 
halo i s not spherically symmc t ric( Aarscth fc Binncy 



1978, Aguilar fc Merritt 199C , Binney 1994) so that 



we study here axially symmetric halo models and cal- 
culate m macho . We consider two models; 
1) Exponential Disk Model 

The axially symmetric density distribution func- 
tion p{R, Z) in cylindrical coordinates is given by 



p(R,Z) 



t R 
p exp(- — 



Z d 



(11) 



where Rd and Zh are scale heights. 
2) Elliptical Model 
p(R, Z) is given by 

po 



p(R,Z) 



(1 



+f) 



(12) 



where a and c describe the ellipticity of the equiden- 
sity surface. This is an axially symmetric version of 
the a model. 

In figure 3 we show jVf MACHO as a function of Rd 
in exponential disk models for various aspect ratios 
[Z h /Rd) of the equidensity surface. We see M^ ACHO 
is - 1.0 x lO n M TLMC f or 0.5 < Z h /R d < 1.0 and it 
increases with the decrease of Zh/Rd for Zh/Rd < 0.5. 

In figure 4a and 4b we show M MACHO in ellipti- 
cal models as a function of a and the aspect ratio 
c/a for a = 2.5 and a = 6.0, respectively. M^ ACHO 
is 8.9 x 10 10 M Q r^ at a=10kpc and c=6kpc for 
a = 2.5 and is 7.08 x \Q w M q t^ at a=22kpc and 
c=13.2kpc for a = 6.0. For large a , M^ ACHO con- 
verges, similarly to the a models. 

4. Discussions 

A deep north Galactic pole pro per motion survey( 



Majewski, Munn & Hawley 1996) suggests that the 
halo is not dynamically mixed but contains a signifi- 
cant fraction of stars with membership in correlated 
stellar streams. If MACHOs are also dynamically 
unmixed, it is possible that the density distribution 
function is neither spherically nor axially symmetric 
but completely inhomogeneous. In such a case what 
we can say from the microlensing events toward LMC 
is the minimal column density S^ ACHO of MACHOs. 

MACHO 



Since in Equation (1), x(l 
is given by 



- x/D s ) < D s /4, E 



iMACHO 



25M pc" 



-2 LMC 
T -6.7 ■ 



(13) 



Similar to Equation (7) , the linear size of the clump of 
MACHOs should be larger than 174pc(a;/kpc) where 
x is the distance to the clump of MACHOs . For 
x = Z) s /2,M^ ACHO is 1.5 x 1O 9 M . If this is the 
case, the optical depth toward the Small Magellanic 
Cloud will be quite different and the inhomogeneity of 
the density distribution of MACHOs can be checked. 

In conclusion, it is shown that the total mass of 
MACHOs becomes smaller than 5 x 1O 1O M , i.e. ~ 
10% of the mass of the halo inside the LMC, either 
if the density distribution of MACHOs is unusual or 
t lmc ^ 2 x 1Q -7 

We like to thank Hayward for checking English. 
This work was supported by Grant-in- Aid of Scientific 
Research of the Ministry of Education, Culture, and 
Sports, No.07640399. 

Appendix 

The proof that M MACHO is minimized if the 
spherically symmetric density distribution func- 
tion p(r) is a constant. 

We first fix the density at r — r c . Since the matter 
inside r c does not contribute to r, M MACHO is mini- 



3 



mized if p(r) is constant for r < r c . For r > r c from 
Equation (4) it is easy to observe that dr/dm is a 
decreasing function of r irrespective of p(r). Now for 
fixed m (i.e. the same mass but different r depend- 



ing! on trie density distribution tunctionj, r is smallest 



and dr/dm is largest if p(r) = p(r c ) =constant. This 
means that r is largest for constant density distribu- 
tion for fixed total mass M MACHO . Inversely for fixed 
r, M MACHO is minimized for the constant density dis- 
tribution. By varying the value p(r c ) , M^^ CHO is 
obtained. 
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Fig. 1— The total mass of MACHOs for the poly- 
tropic model in the unit of Mq . N is the polytropic 
index and R p is the radius of the density distribution 
of MACHOs. 

Fig. 2. — The total mass of MACHOs for the a model 
in the unit of Mq . The density distribution is given 



by p(r) = 




Fig. 3.— The total mass of MACHOs for the expo- 
nential disk model in the unit of Mq. The density 
distribution is given by p(R, Z) = po exp(— — ^-). 

Fig. 4.— The total mass of MACHOs for the elliptic 
model in the unit of Mq . The density distribution 
is given by p(R, Z) = - — -f" 2 — . a) for a — 2.5 and 

b) for a — 6.0 
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